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ABSTRACT 

Multidimensional cosmologies allow for variations of fundamental physical constants 
over the course of cosmological evolution, and different versions of the theories predict 
different time dependences. In particular, such variations could manifest themselves as 
changes of the proton-to-electron mass ratio fi = m p /m e over the period of ~ 10 10 yr 
since the moment of formation of high-redshift QSO spectra. Here we analyze a new, 
high-resolution spectrum of the z = 2.81080 molecular hydrogen absorption system 
toward the QSO PKS 0528-250 to derive a new observational constraint to the time- 
averaged variation rate of the proton-to-electron mass ratio. We find fi\ < 1.5x10" 14 
yr -1 , which is much tighter than previously measured limits. 



Subject headings: cosmology: observations - quasars: absorption lines - quasars: indi- 
vidual (PKS 0528-250) - atomic data 
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1. INTRODUCTION 

The possibility of the variability of fundamental 
physical constants was first put forward by Dirac 
(1937) in the course of his discussion with Milne 
(1937). Later it was considered by Teller (1948), 
Gamow (1967), Dyson (1972) and other physicists. 
Interest in the problem increased greatly during the 
last decade, due to new developments in the Kaluza- 
Klcin and supergravity models of unification of all 
the physical interactions. Chodos & Detweiler (1980), 
Freund (1982), Marciano (1984), and Maeda (1988) 
discussed possibilities of including these multidimen- 
sional theories into the cosmological scenario of the 
expanding Universe and found that the low-energy 
limits to the fundamental constants might vary over 
the cosmological time. Variations of the coupling con- 
stants of strong and electroweak interactions might 
then cause the masses of elementary particles to 
change. Note that an increase of the proton mass by 
0.08% would lead to transformation of protons into 
neutrons (by electron capture), resulting in destruc- 
tion of atoms in the Universe. As demonstrated by 
Kolb, Perry, & Walker (1987) and Barrow (1987), ob- 
servational bounds on the time evolution of extra spa- 
tial dimensions in the Kaluza-Klein and superstring 
theories can be obtained from limits on possible vari- 
ations of the coupling constants. Damour & Polyakov 
(1994) have developed a modern version of the string 
theory which assumes cosmological variations of the 
coupling constants and hadron-to-electron mass ra- 
tios. Therefore the parameters of the theory can be 
restricted by testing cosmological changes of these ra- 
tios. 

The present value of the proton-to-electron mass 
ratio is fi = 1836.1526645 (57) (CODATA, 1997). Ob- 
viously, any significant variation of this parameter 
over a small time interval is excluded, but such vari- 
ation over the cosmological time ~ 1.5 x 10 10 yr re- 
mains a possibility. This possibility can be checked 
by analyzing spectra of high-redshift QSOs. 

The first analysis of this kind has been performed 
by Pagel (1977), who obtained a restriction \fi/fj,\ < 
5 x 10~ n yr -1 on the variation rate of n by com- 
parison of wavelengths of H I and heavy- ion absorp- 
tion lines, as proposed by Thompson (1975). This 
technique, however, could not provide a fully conclu- 
sive result, since the heavy elements and hydrogen 
usually belong to different interstellar clouds, moving 
with different radial velocities. In this paper we use 



another technique, based on an analysis of H2 absorp- 
tion lines only. 

One object suitable for such analysis is the z = 
2.811 absorption system toward PKS 0528-250, in 
which Levshakov & Varshalovich (1985) identified 
molecular hydrogen absorption lines based on a spec- 
trum obtained by Morton et al. (1980). Foltz, Chaf- 
fee, & Black (1988) have presented a limit to pos- 
sible variation of [i based on their observations of 
PKS 0528-250. Their analysis did not, however, 
take into account wavelength-to-mass sensitivity co- 
efficients, hence their result appeared to be not well 
grounded. Subsequently the spectrum of Foltz, Chaf- 
fee, & Black (1988) was reappraised by Varshalovich 
& Levshakov (1993), who obtained \Afi/fi\ < 0.005 
at the redshift z = 2.811, and by Varshalovich & 
Potekhin (1995), who obtained \Afi/fi\ < 0.002 at the 
2a significance level. (Here Afi/fi is the fractional 
variation of fi.) More recently, Cowie & Songaila 
(1995) used a new spectrum of PKS 0528-250 ob- 
tained with the Keck telescope to arrive at the 95% 
confidence interval —5.5 x 10~ 4 < A/i//i < 7 x 10~ 4 . 

Here we present a profile fitting analysis of a 
new, high-resolution spectrum of PKS 0528—250, ob- 
tained in November 1991 with the Cerro-Tololo Inter- 
American Observatory (CTIO) 4 m telescope. We 
have calculated the wavelength-to-mass sensitivity co- 
efficients for a larger number of spectral lines and em- 
ployed them in the analysis, which yields the strongest 
observational constraint yet to possible /x variation 
over the cosmological time scale (eq. fl below). 

2. OBSERVATIONS 

Observations were obtained with the CTIO 4 m 
telescope in a series of exposures of typically 2700 
s duration totaling 33750 s duration. The CTIO 
Echelle Spectrograph with the Air Schmidt camera 
and Reticon CCD was used at the Cassegrain focus 
to obtain complete spectral coverage over the wave- 
length range A w 3465 — 4905 A. Observations of stan- 
dard stars and of a Th-Ar comparison arc lamp were 
obtained at intervals throughout each night, and ob- 
servations of a quartz lamp were obtained at the be- 
ginning or end of each night. For all observations, the 
slit was aligned to the parallactic angle. 

Data reduction was performed following proce- 
dures similar to those described previously by Lanzetta 
et al. (1991). One-dimensional spectra were extracted 
from the two-dimensional images using an optimal 
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extraction technique, and individual one-dimensional 
spectra were coadded using an optimal coaddition 
technique. Wavelength calibrations were determined 
from two-dimensional polynomial fits to spectral lines 
obtained in the Th-Ar comparison arc lamp observa- 
tions. Continua were fitted to the one-dimensional 
spectra using an iterative spline fitting technique. 

The spectral resolution was measured from spec- 
tral lines obtained in the Th-Ar comparison arc lamp 
exposures. This is appropriate because for all ob- 
servations the seeing profile was larger than the slit 
width. The spectral resolution was found to be 
FWHM « 21 - 24 km s _1 in the spectral intervals 
used for the analysis. 

Figure |l| presents parts of the spectrum, with the fit 
superposed on the data, for several spectral intervals 
in which the H2 absorption lines have been analyzed 
(for more detail, see Circovic et al., 1998). 

3. SENSITIVITY COEFFICIENTS 

The possibility of distinguishing between the cos- 
mological redshift of the spectrum and wavelength 
shifts due to a variation of fj, arises from the fact that 
the electronic, vibrational, and rotational energies of 
a molecule each undergo a different dependence on 
the reduced mass of the molecule. Hence comparison 
of the wavelengths of various electronic-vibrational- 
rotational molecular absorption lines observed in the 
spectrum of a high-redshift quasar with correspond- 
ing molecular lines observed in the laboratory may 
reveal or limit the variation of /i with time. 

If the value of \x at the early epoch z of the ab- 
sorption spectrum formation were different from the 
contemporary one, then the ratio 

would deviate from unity. Here 

Ki = dlnAj/dln^ (2) 

is the coefficient which determines the sensitivity of 
the wavelength Xi of ith spectral line with respect to 
the variation of the mass ratio /x. 

These coefficients were calculated previously by 
Varshalovich & Levshakov (1993) from the spectro- 
scopic constants of the H2 molecule, using the Born- 
Oppenheimer approximation. Later Varshalovich & 
Potekhin (1995) calculated Ki in another way, by 



comparison of the H2 laboratory wavelengths with the 
corresponding wavelengths for D2 and T2 molecules, 
which simulate just the mass variation of the study, 
and also for HD molecules. Varshalovich & Potekhin 
(1995) also removed some inaccuracies from the table 
of Varshalovich & Levshakov (1993). The two ways 
of performing the calculation yielded very similar Ki 
values, which argues that both are correct. 

For each electron- vibration-rotational band, a wave- 
length of a transition between two states with the vi- 
brational and rotational quantum numbers v, J and 
v' , J' can be presented as 



A = [v^j, - V l v „j„] 



(3) 



where v is the level energy in cm -1 , and the super- 
scripts u and I stand for the upper and the lower level, 
respectively. For each of them 
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(j(j + i)y 



(4) 



We consider the Lyman bands (transitions X x £+ — > 
B 1 E+) and the Werner bands (X 1 !]^ -> 0*11+) of 
the molecular H2 spectrum. The parameters Y mn for 
the three corresponding states are taken from Huber 
& Herzberg (1979). The coefficient Y 00 is redefined so 
that the energy of each rotational- vibrational band is 
counted from the ground-state energy. For the state 
n u , the factor J( J + 1) in the terms with n = 1 of 
equation (0) has been replaced by (J( J + 1) — 1), in 
order to take into account the projection (A 2 = 1) of 
the electron orbital moment on the molecular axis. 

From the Born-Oppenheimer approximation we 
conclude that the coefficients Y mn are proportional 
to /i~™~ m / 2 . Then the sensitivity coefficients Ki are 
easily obtained from equations (j^)-(^): 



K v , 



J>-v"J 



I — A,/ 



J'-V"J' 



K v"J" 



(5) 



k v J=J2y™n(v + l) (J(J + 1))", (6) 



where 



and the coefficients y mn are given in Table ^ (in 
cm -1 ). For the state n„, the factor J (J + 1) in the 
terms with n = 1 has been replaced by ( J( J+ 1) — 1), 
as well as in equation (El). 
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Fig. 1. — Selected parts of the spectrum PKS 0528—250 with absorption lines of H2 molecules at the redshift 
z = 2.81080. The fit is superposed on the data. The most distinct absorption lines of the Lyman band are labeled 
on the plot. The thick bars along the horizontal axes mark the spectral intervals used in the fit (Sect. 4.1), and the 
asterisks mark the positions of individual lines listed in Table || and used in the independent analysis in Section 4.2 
(note that there are other spectral intervals and lines included in the analyses but belonging to Lyman and Werner 
branches not shown in the figure). 
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4. RESULTS OF ANALYSIS 
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4.1. Synthetic spectrum analysis 
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We have applied to the spectrum a routine de- 
scribed previously by Lanzetta & Bowen (1992). This 
routine performs a comparison of the synthetic and 
observed spectrum and finds an optimal solution to 
a parameterized model of a set of absorption profiles, 
simultaneously taking account of all observed spec- 
tral regions and transitions. Parameter estimates are 
determined by minimizing x 2 , and parameter uncer- 
tainties and correlations are determined by calculat- 
ing the parameter covariance matrix at the resulting 
minimum. 

A total of 59 H2 transitions are incorporated into 
the x 2 fit, and the absorption lines corresponding to 
these transitions occur across the linear, saturated, 
and damped parts of the curve of growth. The red- 
shift, Doppler parameter, and column densities of the 
H2 rotational levels J" = through J" = 7 were 
adopted as free parameters. Wavelengths, oscillator 
strengths, and damping coefficients of the H2 transi- 
tions were taken from Morton & Dinerstein (1976). 
According to Circovic et al. (1998), the total H2 col- 
umn density is logiV(H 2 ) = 18.45 ± 0.02 and the 
Doppler parameter is b — 3.23 ± 0.11 at the redshift 

z = 2.8107998(24). (7) 

Full details of the reduction and analysis of the spec- 
trum are described in a companion paper (Circovic 
et al., 1998), including the list of all spectral inter- 
vals and transitions used in the fit. These spectral 
intervals (shown in Fig. [l] by horizontal bars) were 
chosen to embrace anticipated positions of distinct 
and presumably unblended H2 lines. We emphasize 
that, although the choice of the window function is 
somewhat arbitrary, it should not entail systematic 
shifts of the parameter estimates. Note that there are 
H2 lines present in the spectrum but not used in the 
fit. Some of them are seen in Figure [l] (for example, 
L 1-0 P(2), L 7-0 P(l) and R(2), and others). The 
wavelengths and strengths of these lines are perfectly 
reproduced by the synthetic spectrum. Furthermore, 
within the errors, none of the model lines drops below 
the observed spectrum. This remarkable agreement 
between the measured and model synthetic spectrum 
confirms the reliability of the derived parameters. 

A limit to the variation of the proton-to-electron 
inertial mass ratio was obtained by repeating the x 2 
synthetic spectrum fitting analysis with an additional 
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Fig. 2. — Best fit (with respect to all other parameters 
of the problem) dependence of \ 2 on A/i//x. There 
are 1367 degrees of freedom in the \ 2 fitting analysis. 

free parameter A/1//1. The dependence of x 2 on this 
parameter is shown in Figure 0. The resulting param- 
eter estimate and la uncertainty is 

A/i//i = (8.3tjj;j») x 10- 5 . (8) 

This result^ indicates a value of A/i//x that differs 
from at the 1.6cr level. The 2tr confidence interval 
to Afi/n is 

- 1.7 x 10~ 5 < Afi/n < 2 x 10" 4 . (9) 

Assuming that the age of the Universe is ~ 15 Gyr 
the redshift z = 2.81080 corresponds to the elapsed 
time of 13 Gyr (in the standard cosmological model 
with Q = 1). Therefore we arrive at the restriction 

|/VmI < 1-5 x 10~ 14 yr" 1 (10) 

on the variation rate of //, averaged over 87% of the 
lifetime of the Universe. 

4.2. Profile analysis of separate lines 

We have analyzed the spectrum also by another, 
traditional, technique. The use of the alternative 
technique provides an independent check for the re- 
sults of the above \ 2 analysis and enables a direct 
comparison with the previous results (Foltz, Chaffee, 
& Black, 1988; Varshalovich & Levshakov, 1993; Var- 
shalovich & Potekhin, 1995; Cowie & Songaila, 1995). 



The estimate (H) has been presented at the 17th Texas Symp. 
on Relativistic Xstrophys. (Varshalovich et al., 1996). 



5 



We have selected spectral lines of the Lyman and 
Werner bands which can be unambiguously fitted by 
a single Gaussian profile and a few lines whose de- 
composition in two contours is quite certain (such as 
L 1-0 P(l) and R(2) at A 2 = 4170A and L 7-0 P(l) 
and W 0-0 P(3) at X z = 3866 A). Since there are 
overlapping diffraction orders, we have selected them 
so as to work far from the order edges, and there- 
fore the resolution in the analyzed regions was rela- 
tively high (R > 10 000) and uniform. The analyzed 
50 transitions are listed in the first column of Ta- 
ble H and marked in Figure [l] by asterisks. Only 26 
of these 50 lines have been included in the analysis 
of the synthetic spectrum described in Section 4.1, so 
that the total number of H2 wavelengths analyzed by 
both techniques amounts to 83. Thus in this section 
we use not only an independent technique but also an 
independent choice of the spectral regions. 

The rest-frame wavelengths Ao adopted from Ab- 
grall et al. (1993) are given in the second column of 
Table The third and fourth columns of Table ^ 
present the optimal vacuum heliocentric position of 
the center of each observed profile (A z ) and the es- 
timated standard deviation (a\). The values of A 2 
and <j\ have been provided simultaneously by the 
standard fitting procedure that minimized root-mean- 
square deviations between the fit and data. Sensitiv- 
ity coefficients Ki, calculated according to equations 
(||), (^), are listed in the fifth column. The last col- 
umn presents the redshift corresponding to each A 2 . 
These redshifts, Z{, are shown by crosses in Figure ||. 
In order to test the influence of possible uncertainty 
in Ao , we have repeated the analysis using a set of the 
wavelengths by Morton & Dinerstein (1976); the cor- 
responding redshifts are shown in Figure I with open 
circles. 

In the linear approximation, z(Ki) = z + bKi, 
where b — (1 + z)A/i//i, and z is the cosmological 
redshift of the H2 system. In order the estimates of 
the regression parameters to be statistically indepen- 
dent, it is convenient to write the regression in the 
form 

Zl =z + b{Ki-K), (11) 

where K is the mean value of Ki. 

Given the typical values of a\ ~ 0.02 A and A ~ 
4000 A, one has a typical relative error a z « (a\ /A) (1+ 
z) ~ 2 x 10~ 5 for an individual line. This estimate is 
only an intrinsic statistical error, and it does not in- 
clude an error due to possible unresolved blends. For 
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b = (i+z )A/j./ l j. 




a: + 

b: o 
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-0.01 0.01 0.02 0.03 0.04 0.05 0.06 
sensitivity coefficient K i 

Fig. 3. — Relative deviations of the redshift values, in- 
ferred from an analysis of separate spectral features, 
plotted vs. sensitivity coefficients. The lines repre- 
sent 2a deviations from the slope b of the best lin- 
ear regression. The results based on the rest-frame 
wavelengths by (a) Abgrall et al. (1993) (crosses and 
solid lines) and (b) Morton & Dinerstein (1976) (cir- 
cles and dashes) are shown. The errorbar to the right 
represents the ±2cr limit on Zq. 



this reason, we have not relied on the estimated o\ 
values in our statistical analysis but calculated the la 
uncertainties from the actual scatter of the data (see 
also a discussion in Potekhin & Varshalovich, 1994). 

The estimated mean and slope parameters of the 
linear regression (^TJ), based on Ao by Morton & 
Dinerstein (1976) are z = 2.8107973 (52) and b = 
(-5.85 ± 2.86) x 10~ 4 . Using the data of Abgrall 
et al. (1993), we obtain z = 2.8108028(53) and 
b = (-4.38 ±2.91) x 10~ 4 . The dashed and solid lines 
in Figure ^| correspond to the 2cr-deviations of b for 
the first and second set of Ao, respectively. The er- 
rorbar to the right represents the 2a limit on zq. 

The latter estimate of b translates into 



A/x/M= (-H-5 ±7.6) x 10" 



(12) 



When using the weights oc a x , we obtain a simi- 
lar estimate, Afx/fi = (-10.2 ± 8.1) x 10~ 5 . Since 
the distribution of random errors caused by different 
sources (including possible blends) is not expected to 
be Gaussian, it may worth using robust statistical 
techniques such as the trimmed-mean regression anal- 
ysis (cf. Potekhin & Varshalovich, 1994). We have 
applied the trimmcd-mcan technique of Ruppert & 
Carroll (1980) and found that, at any trimming level 
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up to 12%, the estimate of b is closer to zero but has 
a larger estimated la error, compared with the re- 
sult of the standard least-square analysis reported in 
equation (|l2|). Thus we adopt equation ( |l2| ) as the 
final result of this section. 

The estimate ([l2]) has a larger statistical error com- 
pared with equation (|J). Within 2c, both estimates 
are consistent with the null hypothesis of no variation 
of (U. 

5. CONCLUSIONS 

We have obtained a constraint to the variation 
rate of the proton-to-electron mass ratio fi. Two 
fitting procedures have been used, one of which si- 
multaneously takes into account all observed spectral 
regions and transitions, while the other is applied 
to each spectral feature separately. The two tech- 
niques, applied to two different sets of spectral inter- 
vals, have resulted in similar upper bounds on A/z//z, 
at the level ~ 2 x 10~ 4 . The obtained restriction 
on ( |l0| ) is by an order of magnitude more strin- 
gent than the limit set previously by Varshalovich & 
Potekhin (1995), who used a spectrum with a lower 
spectral resolution. Moreover, it is much more restric- 
tive than the estimate of Cowie & Songaila (1995), 
based on high-resolution Keck telescope observations. 
There are two reasons for the higher accuracy of the 
present estimate. First, our fitting procedure simul- 
taneously takes into account all observed spectral re- 
gions and transitions. This is particularly important 
because many of the transitions are blended, even 
at the spectral resolution of the spectrum used by 
Cowie & Songaila (1995). A separate analysis of spec- 
tral lines leads to larger statistical errors, as we have 
shown explicitly in Section 4.2. Second, we include a 
larger number of transitions between excited states of 
the H2 molecule (83 spectral lines, compared with 19 
lines used by Cowie & Songaila), many of which have 
higher wavelength-to-mass sensitivity coefficients Ki. 
The larger interval of Ki values results in a higher 
sensitivity to possible mass ratio deviations. 

The method used here to determine the variation 
rate of (x could be formally less sensitive than the one 
based on an analysis of relative abundances of chem- 
ical elements produced in the primordial nucleosyn- 
thesis (Kolb et al., 1986). However, the latter method 
is very indirect since it depends on a physical model 
which includes a number of additional assumptions. 
Therefore the present method seems to be more reli- 



able. 

Quite recently, Wiklind & Combes (1997) used a 
similar method (following Varshalovich & Potekhin, 
1996) in order to infer limits on time variability of 
masses of molecules CO, HCN, HNC and the molecu- 
lar ion HCO + from high-resolution radio observations 
of rotational lines in spectra of a few low-redshift 
(z < 1) quasars. The result reported in this paper 
constrains the mass of the H2 molecule, and thus the 
proton mass, at much larger z. These constraints may 
be used for checking the multidimensional cosmologi- 
cal models which predict time-dependences of funda- 
mental physical constants. The described method of 
the calculation of the sensitivity coefficients can also 
be used for analyzing any other high-redshift molecu- 
lar clouds, which may be found in future observations. 

AYP, AVI and DAV acknowledge partial support 
from RBRF grant 96-02-16849 and ISF grant NUO 300. 
KML was supported by NASA grant NAGW-4433 
and by a Career Development Award from the Dud- 
ley Observatory. 

REFERENCES 

Abgrall, H., Roueff, E., Launa, Y. F., Roncin, J. Y., 
Subtil, J. L., 1993, A&AS, 101, 273; ibid., 323 

Barrow, J.D., 1987, Phys. Rev. D, 35, 1805 

Chodos, A., & Detweiler, S., 1980, Phys. Rev. D, 21, 
2167 

Circovic, M., Lanzetta, K. M., Baldwin, J., Williger, 
O, Carswell, R. F. 1998, ApJ (submitted) 

CODATA Bull., Jul. 24, 1997 

Cowie, L. L., & Songaila, A. 1995, ApJ, 453, 596 

Damour, T., & Polyakov, A. M. 1994, Nucl. Phys. B, 
423, 532 

Dirac, P. A. M., 1937, Nature, 139, 323 

Dyson, F. J. 1972, in A. Salam, E. P. Wigner (eds.), 
Aspects of Quantum Theory (Cambridge: Cam- 
bridge Univ. Press), p. 213 

Foltz, C. B., Chaffee, F. H., & Black, J. H. 1988, ApJ, 
324, 267 

Freund, P., 1982, Nucl. Phys., B 209, 146 



Gamow, O, 1967, Phys. Rev. Lett., 19, 759 



7 



Hubcr, K. P., & Herzberg, G. 1979, Molecular Spec- 
tra and Molecular Structure. IV. Constants of Di- 
atomic Molecules (Princeton: Van Nostrand) 

Kolb, E. W., Perry, M. J., Walker, T. P. 1986, Phys. 
Rev. D, 33, 869 

Lanzetta, K. M., & Bowcn, D. V. 1992, ApJ, 391, 48 

Lanzctta, K. M., Wolfe, A. M., Turnshck, D. A., Lu 
L., McMahon, R. G., & Hazard, C. 1991, ApJS, 
77, 1 

Levshakov, S. A., & Varshalovich, D. A. 1985, MN- 
RAS, 212, 517 

Maeda, K.-L, 1988, Mod. Phys. Lett., A 3, 243 

Marciano, W. J., 1984, Phys. Rev. Lett., 52, 489 

Milne, E., 1937, Proc. R. Soc, A 158, 324 

Morton, D. C, & Dinerstein, H. L. 1976, ApJ, 204, 1 

Morton, D. C, Jian-sheng C, Wright A. E., Peterson 
B. A., & Jaunsley, D. L. 1980, MNRAS, 193, 399 

Pagel, B., 1977, MNRAS, 179, 81P 

Potckhin, A.Y., & Varshalovich, D.A. 1994, A&AS, 
104, 89 

Ruppert, D., & Carroll, P. 1980, J. Am. Stat. Ass., 
75, 828 

Teller, E., 1948, Phys. Rev., 73, 801 

Thompson, R., 1975, Astron. Lett., 16, 3 

Varshalovich, D. A., & Levshakov, S. A. 1993, JETP 
Lett., 58, 237 

Varshalovich, D. A., & Potckhin, A. Y. 1995, Space 
Sci. Rev., 74, 259 

Varshalovich, D. A., & Potckhin, A. Y. 1996, Astron- 
omy Lett., 22, 1 

Varshalovich, D. A., Lanzetta, K. M., Baldwin, J. A., 
Potekhin, A. Y., Williger, G. M., & Carswell, R. F. 
1996, in 17th Texas Symp. on Relativistic Astro- 
phys., MPE report 261, ed. W. Vogcs, G. Wieden- 
mann, G. E. Morrill, & J. Triimper (Garching: 
MPE) 265 

Wiklind, T., & Combes, F. 1997, A&A, 328, 48 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v4.0. 



Table 1: Coefficients y mn , in cm 







B 1 ^ 


c 1 ^ 


yio 


2200.607 


679.05 


1221.89 


2/20 


-121.336 


-20.888 


-69.524 


2/30 


1.2194 


1.0794 


1.0968 


2/40 




-0.1196 


-0.0830 


2/50 




0.00540 




2/oi 


60.8530 


20.01541 


31.3629 


2/ii 


-4.0622 


-1.7768 


-2.4971 


2/21 


0.1154 


0.2428 


0.0592 


2/31 


-0.0128 


-0.0293 


-0.00740 


2/41 




0.00138 




2/02 


-0.0942 


-0.03250 


-0.0446 


2/12 


0.00685 


0.005413 


0.00185 


2/22 


-0.0012 


-0.0006867 




2/32 




4.148 x 10~ 5 





2,03 f -38 x 10- 4 
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Tabic 2: H2 lines and sensitivity coefficients 



Line 



Ao (A) A z (A) a x (A) 



K x 



zx 



L 0-0 R(l 
L 0-0 R(0 
L 1-0 P(2 
L 1-0 R(2 
L 1-0 P(l 
L 1-0 R(l 
L 1-0 R(0 
L 2-0 R(3 
L 2-0 P(2 
L 2-0 R(2 
L 2-0 P(l 
L 2-0 R(l 
L 2-0 R(0 
L 3-0 R(3 
L 3-0 P(2 
L 3-0 R(2 
L 3-0 P(l 
L 4-0 R(3 
L 4-0 P(2 
L 4-0 R(2 
L 4-0 P(l 
L 4-0 R(l 
L 4-0 R(0 
L 5-0 R(4 
L 5-0 P(3 
L 6-0 P(3 
L 6-0 R(3 
L 7-0 R(2 
W 0-0 P(3) 
L 7-0 P(l 
L 7-0 R(l 
W 0-0 Q(2) 
W 0-0 Q(l) 
W 0-0 R(2) 
L 9-0 R(2) 
L 9-0 P(l) 
L 9-0 R(l) 
L 9-0 R(0) 
W 1-0 R(3) 
W 1-0 R(2) 
L 10-0 P(2) 
L 11-0 P(3) 
L 11-0 R(2) 
L 12-0 R(3) 
W 2-0 R(3) 
W 2-0 R(2) 
L 13-0 P(3) 
L 13-0 R(2) 
L 15-0 P(3) 
L 15-0 R(2) 



1108.633 
1108.127 
1096.438 
1094.244 
1094.052 
1092.732 
1092.195 
1081.712 
1081.267 
1079.226 
1078.923 
1077.697 
1077.140 
1067.474 
1066.899 
1064.993 
1064.606 
1053.977 
1053.283 
1051.498 
1051.031 
1049.964 
1049.367 
1044.542 
1043.501 
1031.192 
1028.983 
1014.977 
1014.504 
1014.326 
1013.436 
1010.938 
1009.771 
1009.023 
993.547 
992.809 
992.013 
991.376 
987.447 
986.243 
984.863 
978.218 
974.156 
967.675 
966.778 
965.793 
960.450 
956.578 
944.331 
940.623 



4224.779 
4222.830 
4178.284 
4169.945 
4169.226 
4164.157 
4162.108 
4122.218 
4120.463 
4112.779 
4111.648 
4106.879 
4104.751 
4068.050 
4065.712 
4058.479 
4057.029 
4016.490 
4013.838 
4007.062 
4005.259 
4001.183 
3998.954 
3980.460 
3976.557 
3929.652 
3921.287 
3867.848 
3866.085 
3865.381 
3862.010 
3852.498 
3848.034 
3845.218 
3786.139 
3783.448 
3780.378 
3777.961 
3762.962 
3758.449 
3753.172 
3727.786 
3712.273 
3687.606 
3684.203 
3680.493 
3660.059 
3645.243 
3598.670 
3584.501 



0.014 
0.010 
0.012 
0.012 
0.012 
0.008 
0.007 
0.011 
0.015 
0.012 
0.008 
0.016 
0.080 
0.095 
0.117 
0.020 
0.019 
0.012 
0.010 
0.010 
0.008 
0.030 
0.014 
0.031 
0.005 
0.050 
0.010 
0.007 
0.014 
0.009 
0.015 
0.010 
0.011 
0.014 
0.032 
0.023 
0.014 
0.042 
0.019 
0.018 
0.009 
0.011 
0.013 
0.029 
0.013 
0.009 
0.017 
0.010 
0.042 
0.029 



-0.00818 
-0.00772 
-0.00453 
-0.00252 
-0.00234 
-0.00113 
-0.00064 
0.00165 
0.00206 
0.00394 
0.00422 
0.00535 
0.00587 
0.00758 
0.00812 
0.00989 
0.01026 
0.01304 
0.01369 
0.01536 
0.01580 
0.01681 
0.01736 
0.01485 
0.01584 
0.02053 
0.02262 
0.02914 
-0.01045 
0.02976 
0.03062 
-0.00686 
-0.00570 
-0.00503 
0.03647 
0.03719 
0.03796 
0.03858 
0.00439 
0.00562 
0.03854 
0.03896 
0.04295 
0.04386 
0.01324 
0.01456 
0.04574 
0.04963 
0.05430 
0.05816 



2.8107969 
2.8107782 
2.8107765 
2.8108000 
2.8108116 
2.8107761 
2.8107772 
2.8108347 
2.8107800 
2.8108598 
2.8108747 
2.8107884 
2.8107940 
2.8108982 
2.8107678 
2.8107963 
2.8108267 
2.8107983 
2.8107950 
2.8108164 
2.8107905 
2.8108029 
2.8108286 
2.8107082 
2.8107950 
2.8107897 
2.8108264 
2.8107740 
2.8107942 
2.8107576 
2.8108155 
2.8108040 
2.8107949 
2.8108064 
2.8107220 
2.8108365 
2.8107804 
2.8107564 
2.8107874 
2.8108636 
2.8108453 
2.8108005 
2.8107582 
2.8107937 
2.8107977 
2.8108508 
2.8107672 
2.8106998 
2.8108142 
2.8107571 
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